PROCEEDINGS 

0F SCIENCE 



AGATA: Gamma-ray tracking in segmented HPGe 
detectors 



P.-A. Soderstrom , A. Al-Adili, J. Nyberg 

QQ Department of Physics and Astronomy, Uppsala University, SE-75121 Uppsala, Sweden 

E-mail: P-A . Soderstrom@physics .uu.se 

o 

^ F. Recchia 

> 

Q INFN Legnaro, 35020 Legnaro (Padova), Italy 

^ E. Farnea 

(N) INFN Padova, 35122 Padova, Italy 

A. Gadea 

X 

^ IFIC CSIC, Valencia, Spain 

The next generation of radioactive ion beam facilities, which will give experimental access to 
i i many exotic nuclei, are presently being developed. At the same time the next generation of high 

resolution y-ray spectrometers, based on y-ray tracking, for studying the structure of these exotic 
nuclei are being developed. One of the main differences in tracking of y rays versus charged 
particles is that the y rays do not deposit their energy "continuously" in the detector, but in a 
few discrete steps. Also, in the field of nuclear spectroscopy, the location of the source is mostly 
well known while the exact interaction position in the detector is the unknown quantity. This 
makes the challenges of y-ray tracking in germanium somewhat different compared to vertexing in 
00 silicon detectors. In these proceedings we present the methods for determining the 3D interaction 
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positions in the detector and how these are used to reconstruct the y-ray tracks in the AGATA 
detector array. We also present preliminary simulation results of a proposed in-beam method to 
measure the interaction position resolution in the germanium detectors. 
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1. Introduction 

It is sometimes said that we now are on the brink of the fourth revolution within nuclear spec- 
troscopy [1]. The first revolution was the discovery of Nal(Tl) scintillator detectors, by which 
one could start to measure quite accurately the energy and intensity of y-ray transitions in various 
radioactive nuclei. These measurements later radically improved when one started using semicon- 
ductor detectors of germanium instead of Nal(Tl) crystals. The second revolution came with the 
possibilities to use in-beam y-ray spectroscopy together with heavy-ion nuclear reactions. This 
made it possible to study specific levels of the nuclei of interest. The latest revolution is the de- 
velopment of sophisticated high-resolution spectroscopy arrays of high-purity germanium (HPGe), 
crowned by EUROBALL [2] and GAMMASPHERE [3]. These gave access to very weak signals 
from high-spin states, unraveling many new nuclear structure phenomena. And now, just around the 
corner, awaits the second generation facilities for radioactive ion beams (RIB) like NuSTAR/FAIR 
[4, 5] and SIPRAL 2 [6] as the fourth revolution. These facilities will make it possible to study 
very short lived exotic nuclei with extreme values of isospin, located in the terra incognita far from 
the line of /3 stability. Indeed there will be very interesting times for nuclear structure physics in 
the years to come. 

These new exotic nuclei will be produced at the RIB facilities with very low cross sections and 
in a high y-ray background environment, which makes the weak y-ray transitions extremely diffi- 
cult to detect with existing spectrometers. In order to distinguish these rare events, new instruments 
with higher efficiency and resolving power are required. To meet the new requirements 45 institutes 
in 12 European countries (as of fall 2008) are collaborating to build the Advanced Gamma Track- 
ing Array (AGATA) [7, 8]. A similar device, the Gamma Ray Energy Tracking Array (GRETA), is 
also being built in the USA [9]. 

2. AGATA 

The AGATA y-ray spectrometer will be built in stages and moved between different host labo- 
ratories in order to maximally exploit the strengths of the different RIB facilities as they come into 
operation. A first version of AGATA, the AGATA Demonstrator consisting of 1/12 of the full array, 
will begin operation at Laboratori Nazionali di Legnaro (LNL) in Italy in the beginning of 2009. 
The final version of AGATA will consist of a HPGe shell with an inner radius of 23 cm, made 
of 180 crystals assembled into 60 triple cluster detectors. The total amount of HPGe in AGATA 
will be 362 kg. Each HPGe crystal has a hexaconical shape and is of the closed-end coaxial type, 
with the electric field along the radial direction. The crystals are electrically segmented into 36 
individual segments, six azimuthal and six axial, plus the central contact (core). The total HPGe 
solid angle coverage will be 82 % of An. The full-energy peak efficiency of the 4n array will be 
~ 45 % for y rays of multiplicity My = 1 and energy 1 MeV, a factor of about four improvement 
over existing arrays, and ~ 25 % for My = 30, which will give several orders of magnitude increase 
of the resolving power compared to existing arrays [7] . 

This gain in detection power is due to novel techniques to reconstruct the tracks of the y rays 
that scatter within the array. In order to do this, a three dimensional position resolution of < 5 mm 
of the interaction points (IPs) within the crystal is required. Such a position resolution conesponds 
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Figure 1: Signals from the core, the segment with the primary hit (Seg 4), and from the mirror charges for 
a y-ray interaction in a six fold segmented HPGe detector. Figure from [7]. 

to an angular resolution of about 1 degree. This is achieved using digital electronics to sample the 
pulse shapes and advanced pulse-processing techniques to locate the IPs. 

2.1 Pulse-shape analysis 

In order to determine the interaction position, pulse-shape information from the segment which 
was hit, the mirror charges in the neighboring segments, and the core signal is used, as illustrated in 
fig. 1. The method to use pulse-shape analysis (PSA) to determine interaction position is, however, 
not unique for AGATA and GRETA. This is for example done also in MINIBALL, where the 
azimuthal position is obtained by comparing the amplitudes of the mirror charges [10]. The radial 
position is obtained from the shape of the signal of the hit segment. This is good enough to correct 
for Doppler effects, but is does not give enough resolution to do y-ray tracking. The PSA algorithm 
that will be implemented in AGATA is the grid search algorithm [11], where the sampled pulse 
shapes are compared to a database of pulse shapes for different interaction positions. 

One of the main problem of PSA algorithms is to obtain the database of pulse shapes. When 
delivered, the crystals are scanned using radioactive sources [12]. This process gives a very precise 
database of pulse shapes but it is unfortunately a very slow process. To complement the scanned 
pulse shapes, codes for calculating them are being developed [13]. These calculations are much 
faster, but the main challenge is to obtain realistic pulse shapes, while dealing with large uncertain- 
ties in the impurity concentration and in the modeling of charge mobility [14]. 

3. Tracking of 7 rays 

Once the interaction points and the corresponding energies are identified, the different y rays 
should be disentangled from the "world map", see fig. 2. Since the primary interaction mechanism 
at these energies is a series of Compton scatterings followed by a photoelectric effect, the tracking 
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Figure 2: Simulated interaction points of 30 y rays of energy Ey = 1 .33 MeV in the (9,<j) sin 9) plane of an 
ideal germanium shell with an inner radius of 15 cm and an outer radius of 24 cm. Circles are correctly and 
squares incorrectly identified clusters. Figure from [7]. 

algorithm is based on the Compton scattering formula 



In order to detect higher and lower energy /-rays efficiently the algorithm should also take into 
account single photoelectric absorption as well as electron-positron pair-production. It is also im- 
portant that the algorithm can discriminate between y rays and other types of interactions, like 
neutrons [15, 16]. Two main tracking algorithms to be implemented in parallel has been developed 
and compared carefully [17]. 

3.1 Clusterisation or Forward tracking 

One of the algorithms to be used is the clusterisation algorithm, developed by the group at 
INFN Padova [18]. In this algorithm the interactions caused by a certain 7 ray are assumed to 
be well localized, which means that the IPs from the same y ray should make up clusters in the 
(0,0sin0) plane, see fig. 2. After these clusters have been identified, the source location is as- 
sumed as the zeroth IP. A first and a second IP is chosen and the y-ray energy after scattering is 
determined from the measured energy depositions in the cluster. This energy is compared to the 
energy calculated from the Compton scattering formula in eq. (3.1) and a figure of merit based on 
the agreement between the two energies is calculated and if necessary the clusters are rearranged. 
All possible permutations are evaluated and the one with the best figure of merit is chosen. This 
procedure is repeated with the identified first IP as a starting point until all IPs in a cluster have 
been assigned to the track. 




1 + Ar(l + COS0) 



(3.1) 
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3.2 Backtracking 

The backtracking algorithm was developed by the group at KTH, Stockholm [19, 20]. This 
algorithm uses the information that the final IP most probable has an energy deposition between 
100 to 250 keV. Starting from this assumed final IP, other IPs are searched for within a distance 
based on the interaction length in germanium for y rays of that energy. This procedure is repeated 
until the track is terminated by the source location. It is then repeated until no more suitable IPs 
are available. 



4. In-beam measurements of position resolution 

As mentioned in sec. 2, the position resolution within the crystal is of fundamental importance 
in order to do y-ray tracking. Thus this is an important parameter to measure experimentally, which 
is planned to be done during the commissioning phase at LNL. One way to do this directly is to 
use a collimated beam of y rays in order to get a well defined first IP. But this requires a very 
narrow beam in order not to be limited by the collimation. The consequence of this would be a 
very time consuming measurement in order to get good statistics. Instead indirect measurements 
based on imaging techniques [21] or in-beam measurements of Doppler correction capabilities are 
being examined. 

The contribution to the full width at half maximum (FWHM) of the full-energy peak in a y-ray 
spectrum (see fig. 3) can be divided into four parts 

W 2 = W? + W* + Wl + Wi y , (4. 1) 

where Wj is the intrinsic resolution of the detector, W Vl and Wq y are the contributions from the uncer- 
tainties in the velocity and angle of the y-ray emitting nucleus (recoil), and Wg is the uncertainty 
in the emission angle of the yray relative to the beam axis. The last three terms in eq. (4.1) are due 
to effects from the Doppler shift, 



1 -j3cos0 _ 1 -/3cos0 7 



E' y = E y — ps Ey — " n J . (4.2) 



where j8 = v r /c and 6 is the angle between the recoil and y-ray directions. The last expression 
is valid if the recoil direction is close to the beam direction (0 degrees). Since the y-ray emission 
angle is determined by the position of the first IP in the detector, there is a direct correspondence 
between We and the position resolution, p, while W\, W Vl and Wg r are independent of p. Since Wg y 
depends on the Doppler correction of the y-ray energy from Ey to Ey as 



dE'y p sin 9y 



W 6y = -^A6 r = Ey ' = A0y, (4.3) 



with Ady being the angular resolution, one would like to have a setup with the detector places as 
close as possible to the target, as high j8 as possible, and the detector placed at an angle close to 90 
degrees. See fig. 4 for an example of a simulation, using the AGATA Geant4 code [22], showing 
how W depends on p for two different nuclear reactions. 
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9Be(80Se,3n)86Sr @ 40 mm, Smearing 5 mm 
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Figure 3: Photo peak from a simulation of the reaction 9 Be( 80 Se,3n) 86 Sr using the clusterisation algorithm 
for tracking, a position resolution of 5 mm and a distance between the detector and the source of 40 mm. 
HereW r 2 = w2+W Q 2. 




Figure 4: Simulations of the total width, W, of the full-energy peak in an AGATA triple cluster detector 
placed at a distance of 15 cm from the target to the front of the detector and at 90 degrees relative to the 
beam as a function of the position resolution. In the left panel the reaction d( 51 V,n) 52 Cr and in the right 
panel the reaction d( 37 Cl,n) 38 Ar is shown. Figure from [23]. 

There is some previous experience using this method to estimate the position resolution of a 
detector. Kroll et al. have done a similar measurement on the MARS detector [24, 25], an Italian 
pre-AGATA detector. The GRETA collaboration has also made measurements of the position res- 
olution of their detectors, using the reaction 12 C( 82 Se,4n) 90 Zr at 385 MeV. They obtained a p of 
4.7 mm [26]. Two experiments have also been made with AGATA detectors, one using a single 
AGATA crystal [27] and another using a prototype triple cluster detector [21, 28]. Both of these 
experiments were performed at the tandem accelerator laboratory at University of Cologne. 

In the experiment with a single AGATA crystal, the fusion-evaporation reaction d( 37 Cl,n) 38 Ar 
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Figure 5: Simulation of the angular and velocity distribution of the recoils in the reaction 9 Be( 80 Se, 3n) 86 Sr. 

at 70 MeV, and in the experiment with a triple cluster detector, the fusion-evaporation reaction 
d( 48 Ti,p) 49 Ti at 100 MeV was used. The results gave a p of 4.3 mm FWHM at a y-ray energy 
of 2168 keV and a p of 4.4 mm FWHM at a /-ray energy of 1382 keV in the two experiments, 
respectively. The use of protons as the evaporated particle complicated the setup and the analysis, 
resulting ultimately very time consuming, time that will not be available during the commissioning 
phase. Furthermore, the main uncertainty in the estimation process was coming from the modeling 
of the reaction by the Monte Carlo simulation used to translate the FWHM of the full-energy peak 
to position resolution in the detector. 

The lessons learned from this experiment is that for commissioning, the experimental setup 
must be simplified even further in order to decrease the analysis time. One also needs to find a 
Monte Carlo independent way to translate a measured FWHM of the full-energy peak to position 
resolution. The way to deal with the first issue is to choose another reaction, that evaporates neu- 
trons instead of protons, since these do not have to cross any Coulomb barrier and thus will have 
lower average energy and the angular deviation of the recoils from degrees with respect to the 
beam will be as small as possible. A possible reaction with such properties is 9 Be( 80 Se,3n) 86 Sr. 
Results of a simulation of this reaction are shown in figs. 3 and 5. As seen, the recoils are well 
localized in the ( r , jS ) plane. 

A new strategy has been proposed in order to obtain a Monte Carlo independent estimation of 
the position resolution [21]. The idea is that instead of taking a measurement at one distance and 
compare it to a Monte Carlo simulation, measurements of the total widths W c i ose and Wf ar are made 
at two distances <i c iose and <if ar . The position resolution can then be estimated through 

^ (^-^(jjL -£)-', (4.4) 

with k being a constant independent of W and d. Such a setup has been simulated for different 
reactions of which preliminary results are shown in fig. 6. The deviation of the estimated values 
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Figure 6: Position resolution as a function of smearing in the tracking program for the reactions 
9 Be( 80 Se,3n) 86 Sr, 2 H( 81 Br,n) 82 Kr and 9 Be( 136 Xe,3n) 142 Ce for the two distances d c i ose = 40 mm and 
dfar =140 mm. 

from the simulated ones is due to the fact that AGATA Demonsuator placed close to the target is 
covering a wide angle around 90 degrees. Another contribution to this deviation is the uncertainty 
in the distance to the front face of the detector, due to its curvature radius, when the detector is 
placed closer to the target than the nominal distance of 235 mm. Future investigations are foreseen 
to refine the estimation procedure. For more details, see ref . [29] . 

5. Summary 

In order to study weak transitions in nuclei with extreme values on spin, isospin or temperature, 
new radioactive ion beam facilities are now being built. These new facilities will be complemented 
by detector arrays with a great increase in resolving power. The novel technique implemented 
in these arrays are that of y-ray tracking using pulse-shape analysis with digital electronics. The 
efficiency of the tracking depends, however, critically on the obtainable position resolution of the 
detector. Some preliminary results on simulations to determine the position resolution of the the 
AGATA Demonstrator y-ray tracking array have been presented. 
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